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ABSTRACT 

We use the clustering of galaxies around distant AGN to show with ~ 90% confidence that fainter AGN 
are longer lived. Our argument is simple: since the measured galaxy-AGN cross-correlation length 
tq ~ hh^ 1 Mpc does not vary significantly over a 10 magnitude range in AGN optical luminosity, faint 
and bright AGN must reside in dark matter halos with similar masses. The halos that host bright 
and faint AGN therefore have similar intrinsic abundances, and the large observed variation in AGN 
number density with luminosity reflects a change in duty cycle. 

Subject headings: galaxies: high-redshift — cosmology: large-scale structure of the universe — quasars: 
general 
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1. INTRODUCTION 

In the famous paper that postulated a link between 
quasars and accreting black holes, Lynden-Bell (1969) 
remarked that the black holes created by quasar accre- 
tion would be gigantic and common, with masses around 
10 8 M Q and a space density similar to that of local galax- 
ies. It was a prescient comment, but Soltan's (1982) re- 
finement of his calculation drew attention to the impor- 
tance of the assumed quasar lifetime. The total accretion 
was sufficient to place a 1O 6 M black hole inside every 
galaxy brighter than M31, Soltan showed, but the ac- 
creted mass might equally well be distributed among a 
smaller number of heavier black holes or a larger num- 
ber of lighter ones. The length tQ of the quasars' lives 
would determine which was the case. Although the un- 
derstanding of black hole formation has advanced enor- 
mously since that time, <q remains a key parameter in 
theoretical models. Our ignorance of it is arguably the 
largest source of uncertainty in the accretion histories of 
supermassive black holes. 

This paper is concerned not with the value of the 
quasar lifetime itself, but rather with the idea that there 
is a single lifetime for accretion onto active galactic nu- 
clei (AGN). It is obviously an oversimplification. The 
duration of a luminous accretion episode is presumably 
affected by the mass of the central black hole, the size 
of the gas supply, the nature of the event that funnels 
gas towards the black hole, the strength and duration of 
dust obscuration, and so on. Our aim is to measure the 
extent to which this produces a systematic dependence 
of the lifetime on the luminosity of the AGN. 

It is easy to convince oneself that such a dependence 
might exist. The extreme accretion associated with the 

1 Based, in part, on data obtained at the W.M. Keck Observa- 
tory, which is operated as a scientific partnership between the Cal- 
ifornia Institute of Technology, the University of California, and 
NASA, and was made possible by the generous financial support 
of the W.M. Keck Foundation. 
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most luminous QSOs is rare and must have a small duty 
cycle (e.g., Martini 2004), while low- level accretion has a 
high enough duty cycle to be observed in approximately 
half of all nearby galaxies (e.g., Ho 2004). As far as we 
know, however, no-one has previously attempted a di- 
rect measurement of the dependence of AGN lifetime on 
luminosity (cf. Merloni 2004, Hopkins et al. 2005). Al- 
though it may seem perverse to try to look for systematic 
differences in the accretion lifetime when the lifetime is 
still uncertain by two orders of magnitude (e.g., Martini 
2004) , in fact (as we show in § El changes in the lifetime 
are much easier to measure than the value of the lifetime 
itself. 

Our approach exploits the well known fact that the 
duty cycle of a population of objects can be inferred from 
its number density and clustering strength (e.g., Adel- 
berger et al. 1998). The reason is simple. In universes 
with hierarchical structure formation, the rarest and 
most massive virialized halos cluster the most strongly 
(e.g., Kaiser 1984), and so the mass and number density 
of the sub-population of halos that contain the objects 
can be deduced from the strength of the objects' clus- 
tering. The duty cycle is equal to the objects' observed 
number density divided by the number density of halos 
that can host them. If clustering measurements indicate 
that AGN reside in halos of mass 10 12 M Q , for example, 
but the number density of AGN is only 1% of the num- 
ber density of halos with M = 10 12 M Q , the duty cycle 
is evidently 0.01. 

Martini & Weinberg (2001) and Haiman & Hui (2001) 
were the first to discuss the technique in detail. Our 
treatment is similar to theirs, except in one important 
respect: we infer the duty cycle from the clustering 
of galaxies around AGN, rather than from the clus- 
tering of AGN themselves. As pointed out by Kauff- 
mann & Haehnelt (2002), the high number density of 
galaxies makes the galaxy-AGN cross-correlation length 
much easier to measure than the AGN auto-correlation 
length. A major additional benefit is that any sur- 
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vey deep enough to detect galaxies around bright high- 
redshift QSOs will inevitably detect faint AGN at the 
same redshifts, increasing the sample size and the lumi- 
nosity baseline over which changes in duty cycle can be 
measured. 

2. DATA 
2.1. Galaxies 

The data we analyzed were taken from our color- 
selected surveys of star-forming galaxies with magnitude 
TZab < 25.5 and redshift 1.8 ^ z H 3.5. A more com- 
plete description of the surveys can be found in Steidel 
et al. (2003), Steidel et al. (2004), and Adelberger et 
al. (2005b). We review only the most important aspects 
here. 

The surveys consist of measured redshifts for 1627 
galaxies with redshift z > 1 in 19 fields scattered around 
the sky (tabled . (These totals exclude any survey fields 
with no detected AGN and include only the galaxies with 
the most certain redshifts.) The size of the fields varies 
but is typically ~ 100-200 square arcmin. The coor- 
dinates of some fields were chosen more-or-less at ran- 
dom, but most fields were centered on a bright QSO or 
group of QSOs. Objects were selected for spectroscopy 
if their U n GlZ colors indicated they were likely to lie in 
the targeted range of redshifts. Our decision to obtain a 
spectrum of an object was influenced only by its U n GlZ 
colors, 1Z magnitude, and spatial position; we were more 
likely to observe objects if they had 23 < 1Z < 24.5, if 
they had colors similar to those expected for AGN, or if 
they lay close to a known AGN, and we rarely observed 
objects whose colors did not satisfy the selection criteria 
of Steidel et al. (2003) and Adelberger et al. (2004). 
The overall redshift distribution of the galaxies in these 
fields is shown in figure H Their distribution of absolute 
magnitudes, calculated from observed broadband colors 
for a concordance cosmology with — 0.3, fl\ = 0.7, 
h = 0.7, is shown in figure|21 

2.2. AGN 

Fifty-seven of the 1684 objects in our spectroscopic 
sample have strong emission in both Lyman-a and CIV 
A1549. We classify these objects as AGN for reasons that 
are discussed in Steidel et al. (2002). Although some of 
our faintest AGN might be misclassified as galaxies be- 
cause their CIV lines are too weak for us to detect, the 
lack of CIV emission in the thousand-object composite 
spectrum of Shapley et al. (2003) shows that these mis- 
classified AGN must be rare. 

Our total sample of AGN was increased to 79 by adding 
the previously known AGN that we deliberately included 
in our survey fields. Since CIV was the only line (aside 
from Lyman-a) detected with reasonable significance in 
every AGN spectrum, we based our redshift assignments 
on it. In their analysis of 3814 QSOs from the Sloan 
Digital Sky Survey, Richards et al. (2002) found that 
CIV was blueshifted on average by 824 km s _1 com- 
pared to Mgll, which they assumed was at the QSO's 
systemic redshift. We accordingly assumed that each of 
our AGN's true redshift was 824 km s _1 redder than the 
peak of CIV emission. Since Richards et al. (2002) re- 
port a scatter in the CIV-Mgll velocity offsets of 500 km 
s , we expect that the uncertainty in our QSO redshifts 
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Fig. 1. — Redshift distributions for the galaxies and AGN in 
our sample. Also shown is the number of galaxy- AGN neighbors, 
defined as the number of galaxy- AGN pairs with angular separation 
60" < 9 < 300" [1.2 < H/(/i- 1 comovingMpc) ^ 6.2] and radial 
separation AZ < 30/i _1 comoving Mpc. 
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Fig. 2. — Distribution of absolute AB magnitude at rest-frame 
1350A for the AGN and galaxies in our spectroscopic sample. No 
corrections for incompleteness have been applied, so these do not 
resemble the true distributions for the underlying populations. 



will be approximately 500 km s _1 . Although the way we 
assign redshifts is better suited to our sample's broad- 
lined AGN, any mistakes in the redshifts of narrow-lined 
AGN are unlikely to affect our conclusions: as we will see, 
the typical redshift error would have to be ~ 3000 km 
s -1 (i.e., ~ 30/t" 1 comoving Mpc) to alter our clustering 
measurements significantly. Figure ^ shows the redshift 
distribution for the 79 AGN. Figure |5] shows their distri- 
bution of velocity FWHM and apparent magnitude. 
Strong emission lines prevented us from calculating 
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a Number of (non-active) galaxies with spectroscopic redshift z > 1. 

b Number of AGN with spectroscopic redshift z > 1 and rest-frame 1350A absolute AB magnitude M1350 > —24 
c Number of AGN with spectroscopic redshift 2 > 1 and Afi35o < —24 

d The field is centered on this QSO, but the QSO itself is excluded from our analysis because we lack a good spectrum. 
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Fig. 3. — Overview of the characteristics of the AGN in our sam- 
ple. Upper left: Histogram of CIV line width. The typical uncer- 
tainty ranges from 10—20%, and is dominated by systematics (e.g., 
continuum placement) for the brightest AGN. Upper right: His- 
togram of absoluate AB magnitude at rest-frame 1350A (M1350). 
The uncertainty in the AB magnitude is ^ 0.2 magnitudes for even 
our faintest objects (e.g., Steidel et al. 2003). Lower left: Rela- 
tionship between CIV line width and apparent AB magnitude at 
rest-frame 1350A. Lower right: M1350 against redshift. Recall that 
the selection bias is severe in our AGN sample, since (for exam- 
ple) we deliberately targeted AGN that were bright and had broad 
emission lines. These panels show the characteristics of our sample 
as selected, not of a fair sample of high-redshift AGN. 



the AGNs' AB magnitudes at rest-frame 1350A directly 
from their broadband magnitudes. Instead we scaled 
each AGN's spectrum to match its observed G and 1Z 
magnitudes, measured the flux density near 1350A, then 



converted to absolute magnitude for a cosmology with 
Oil/ = 0.3, Qa — 0.7, h — 0.7. This procedure failed for 
our brightest sources, G ^ 18, which were saturated in 
our images. For these we adopted the magnitude implied 
by their unsealed flux-calibrated spectra. Three of our 
sources were saturated and lacked flux-calibrated spec- 
tra. The magnitudes of these were taken from the Sloan 
Digital Sky Survey archive or from photographic mea- 
surements in the NASA/IPAC Extragalactic Database. 
Figure [21 shows the resulting histogram of AGN absolute 
magnitude. Although unintended, our selection strat- 
egy has given us a sample of AGN with brightnesses dis- 
tributed almost uniformly over a 10-magnitude range. 
Comparison to the galaxies' apparent magnitude distri- 
bution suggests that stellar light may contribute sig- 
nificantly to the measured magnitudes of the faintest 
AGN. We do not correct for this. Doing so would only 
strengthen our conclusions, since the faintest AGN would 
be even fainter than we assume. 

2.3. Simulations 

In a number of places our interpretation of the data re- 
lies on the GIF-LCDM numerical simulation of structure 
formation in a cosmology with Qm = 0.3, Qa — 0.7, h — 
0.7, r = 0.21, and erg = 0.9. This gravity-only simulation 
contained 256 3 particles with mass 1.4 x lO lo /i -1 M in 
a periodic cube of comoving side- length 141.3ft. -1 Mpc, 
used a softening length of 20ft" 1 comoving kpc, and 
was released publicly, along with its halo catalogs, by 
Frenk et al. astro-ph/0007362. Further details can be 
found in Jenkins et al. ("1998) and Kauffmann et al. 
(1999). Although the simulation does not include much 
of the physics associated with galaxy formation, we make 
use only of its predictions for the statistical distribu- 
tion of dark matter on large ( ^, Mpc) scales. Since 
the GIF-LCDM cosmology is consistent with the Wilkin- 
son Microwave-Anisotropy Probe results (Spergel et al. 
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2003), and since modeling the gravitational growth of 
perturbations on large scales is not numerically challeng- 
ing, the large-scale distribution of dark matter in this 
simulation should closely mirror that in the actual uni- 
verse. 

3. METHODS 
3.1. Estimating ro 

We estimated the correlation lengths of the samples 
with two approaches. Both correct for the irregular an- 
gular sampling of our spectroscopy and are unaffected 
by the selection criteria that were used to include AGN 
in our sample. The second approach is also insensitive 
to the criteria that were used to select the galaxies. See 
Adelberger (2005) for a more complete discussion. 

In the first approach, we cycle through the AGN in our 
sample, calculating for each one both the number iV bs(^) 
of galaxies in the AGN's field whose comoving radial sep- 
aration from the AGN, AZ, is less than £ = 30/i _1 Mpc, 
and the number N cxp (£,ro) that would be expected if 
the correlation function had the form £(r) = (r/ro) -1 ' 6 . 
The quantity N cxp (£,r ) is related straightforwardly to 
the integral of the correlation function along the lines of 
sight to galaxies in the field. As shown by Adelberger 
(2005), 

where the sum runs over all galaxies in the AGN's field, 
Zi is the AGN's redshift, Az is the redshift difference 
corresponding to a comoving radial separation of size £, 
Pj(z) is the selection function for the jth galaxy 3 , nor- 
malized so that J Q °° dzPj(z) = 1, and r.y is the distance 
between the AGN and a point at redshift z with the 
galaxy's angular separation 9j. We then sum the val- 
ues N b s (£) and N exp (£,ro) for all our AGN, and take 
as our best-fit correlation length the value of ro that 
makes the total expected neighbor counts equal to the 
total observed. To ensure that our estimate of ro reflects 
the clustering strength on large ( £ Mpc) scales, rather 
than conditions inside the AGNs' halos, we exclude from 
consideration any galaxy-AGN pairs with angular sepa- 
ration 9 < 60" (i.e., X.lhr 1 comoving Mpc at z = 2.5). 
Galaxy-AGN pairs with 6 > 300" are also excluded, since 
the weak clustering signal at the largest angular separa- 
tions can be overwhelmed by low-level systematic errors 
(Adelberger 2005). 

The approach of the preceding paragraph can fail if the 
assumed selection functions Pj are inaccurate. To guard 
against this possibility we also estimate ro by finding the 
value that makes 

X)agn A r Q bs(^) _ Eagn N cxp (£, rp) . . 

Eagn ^obs(2^) ~ Eagn N atp (2t, ro) ' 
Taking the ratio causes most systematic errors to cancel 
(Adelberger 2005). Since it also increases the random 

3 Since the galaxies in our samples were chosen with different 
color-selection criteria, their expected redshift distributions are dif- 
ferent. In this approach, we set Pj to the observed LBG redshift 
distribution if the object was selected with the LBG selection crite- 
ria and to the observed BX redshift distribution if the object was 
selected with the BX criteria. Otherwise the galaxy is ignored. 
(See Adelberger et al. 2004 for a definition of these criteria and 
plots of their redshift distributions.) 



errors, however, we use the equation [5] only to verify 
that systematic errors have not badly compromised the 
estimate of ro from the first approach. 

3.2. Estimating the duty cycle 

As stated in the introduction, our definition of duty 
cycle is the observed number density of AGN divided by 
the number density of halos that can host them. Calcu- 
lating it requires two steps. 

Halo abundance 

We use the GIF-LCDM simulations to estimate the 
halo abundance from ro. For each of the publicly re- 
leased catalogs of halos at redshifts 2 < z < 3 4 , we cal- 
culated the cross-correlation function £mi,m 2 ( r ) of halos 
in two mass ranges, M > M\ and M > M2, for dif- 
ferent choices of Mi and Mi, and estimated the cross- 
correlation length ro by fitting a power-law to £mi,M 2 
at separations 1 < r/(/i -1 Mpc) < 10. After calculat- 
ing the number density of halos with M > M2 in the 
simulations at redshift z, we stored our results as a table 
ro(z, Mi, 712) giving the expected cross-correlation length 
at redshift z between halos with mass threshold Mi and 
halos with number density n^.. If all our observations 
were at redshift zq and we knew the threshold mass M g 
of the galaxies' halos, we could convert any measured 
correlation length ro into a number density n q of AGN 
halos by simply looking up the value of n q that made 
the tabulated ro(zo,M g ,n q ) equal our observed correla- 
tion length. In fact our observations are at a range of 
redshifts and the galaxy mass is not precisely known. 
Figure 0] shows the uncertainty in the relationship be- 
tween ro and n q that results from the range of redshifts 
in our survey and from the la uncertainty in the galaxy 
masses (Adelberger et al. 2005a). For the remainder 
of the paper we will adopt an ro~n q relationship that is 
a least-squares fit to the data in the figure (solid line). 
Although we can offer little justification for this com- 
promise, the exact choice of relationship has almost no 
effect on our conclusions. Any errors in the relationship 
increase or decrease in tandem the implied duty cycles 
for bright and faint AGN; they alter the absolute value 
we infer for the duty cycles but not the relative differ- 
ence between them. (We demonstrate that this is true 
in Figure|B| below.) This is one of the main strengths of 
our approach. It justifies our claim in Sj^that a system- 
atic variation of AGN lifetime with luminosity is easier 
to measure than the absolute value of the lifetime itself. 

AGN abundance 

We adopt a crude approach since small (tens of per- 
cent) errors in the AGN abundance have little effect on 
our conclusions. At the faintest magnitudes we estimate 
the AGN number density by multiplying the galaxy lu- 
minosity function at z = 3 (Adelberger & Steidel 2000) 
by /(M1350), the fraction of sources in our spectroscopic 
sample with absolute magnitude Afi35o that were ob- 
served to be AGN. Note that we are including all AGN 
in this analysis, not merely the broad-lined AGN consid- 
ered by Hunt et al. (2004). Since the faint end of the 
rest-frame UV luminosity distribution of galaxies does 

4 i.e., for the catalogs at z = 2.97, 2.74, 2.52, 2.32, 2.12 
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Fig. 4. — Theoretical relationship between cross-correlation 
length rn and AGN-halo comoving abundance n. Points show the 
GIF-LCDM relationship at two redshifts. The error bars indicate 
the uncertainty in the relationship due to the uncertainty in the 
galaxies' threshold mass. The solid line shows the least-squares 
compromise that we adopt throughout: log(ra/(/i -1 Mpc) 3 ) = 
— 0.83r'o + 1.00. The upper and lower dotted lines show the re- 
lationships that would result if we altered the assumed threshold 
mass by ±la\ Figure Irjl shows that our conclusions would not be 
significantly affected if we adopted these relationships instead. 
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Fig. 5. — Observed number density vs. magnitude for AGN at 
z ~ 2. Squares show the 2dF QSO luminosity function of Croom et 
al. (2004). Circles show our rough estimate of the AGN luminosity 
function at fainter ma gnitud es, calculated from our survey with the 
method described in § 13,2.01 The crude completeness corrections of 
this approach yield a luminosity function adequate only for cases 
like ours where low accuracy is tolerable. The parameters of the 
Schechter-function (solid line; M 4 = —26.2, a = —1.85, = 
4 X 10~ 7 Mpc -3 ) should not be used in other situations. 



not evolve significantly from z — 3 to z — 2 (N.A. Reddy 
et al. 2005, in preparation), this number density should 
be roughly appropriate down to z = 2. At the brightest 
magnitudes we adopt the "2dF" 1.81 < z < 2.10 QSO 
luminosity function of Croom et al. (2004). 5 The AGN 
luminosity distribution is fit tolerably well by a Schechter 
function (figure [SJ) , and we use this fit to estimate the 
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5 We convert the absolute magnitudes M;, 7 reported by Croom et 
(2004) to M1350 by adding 0.46 magnitudes; subtracting 0.07 



magnitudes converts to the AB system, and adding 0.53 undoes 
their /('-correction from observed-frame to rest-frame bj (Cristiani 
k, Vio 1990). 



10° 

O 

|,10- 3 

o 

|-io- 3 

XJ 

io- 4 
io- 5 

-30 -28 -26 -24 -22 -20 
M 

Fig. 6. — Top panel: galaxy-AGN cross-correlation length as a 
function of AGN luminosity M1350. Points with error bars show 
our measurements. The shaded square shows the ±1ct range of the 
galaxy-galaxy correlation length at similar redshifts (Adelberger et 
al. 2005a); its abscissa is arbitrary. Bottom panel: implied duty 
cycle as a function of AGN luminosity. Error bars show the ±1(T 
random uncertainty. The four- and six-pointed stars show how our 
estimated duty cycle would change if we altered the assumed rela- 
tionship between clustering strength and abundance by an amount 
similar to its uncertainty. (They correspond to the upper and lower 
dotted envelopes in Figure [I]) Note that the confidence intervals 
shown in this plot reflect only the constraints from our clustering 
analysis. Other considerations rule out a duty cycle of ^ 1 for the 
faint AGN and < IO" 5 for the bright AGN, however. See §|!>]for 
further discussion. 



number density of AGN in each range of apparent mag- 
nitude. 

4. RESULTS 

The first approach of § 13.11 leads to the estimates 
ro = 4.7, 5.4 comoving Mpc for the galaxy-AGN cross- 
correlation length of AGN with magnitude —30 < 
-^1350 < —25 and —25 < Afi35o < —19, respectively. An 
easy way to estimate the uncertainty is suggested by the 
similarity of the cross-correlation length to the galaxy- 
galaxy correlation length reported by Adelberger et al. 
(2005a): generate many alternate realizations of the data 
by treating randomly chosen galaxies in each field as that 
field's AGN, rather than the true AGN themselves, and 
recalculate ro for each simulated sample. The rms dis- 
persion of ro among these simulated samples should be 
roughly equal to the uncertainty in r , and we adopted it 
for the error bars in the top panel of Figure The true 
uncertainty is likely to be somewhat smaller, since our 
spectroscopic selection strategy gave our AGN more an- 
gular neighbors with measured redshifts than the typical 
galaxy. 

The bottom panel of Figurc[|)]shows the same data, ex- 
cept the cross-correlation length has been converted to a 
duty cycle with the approach of § 13.21 As emphasized in 
that section, uncertainties in the rrj-abundance relation- 
ship mean that the labels on the y axis could be wrong 
by a multiplicative constant, but relative differences in 
duty cycle should be secure. 

To estimate the significance of the apparent difference 
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in duty cycle, we note that our adopted relationship be- 
tween ro and halo number density implies that ro would 
be 2.92ft, -1 comoving Mpc larger for AGN with —30 < 
M 1350 < -25 than AGN with -25 < Mi 350 < -19 un- 
der the null hypothesis that the duty cycle is independent 
of M1350. The observed difference in best- fit correlation 
length, —0.72ft" 1 comoving Mpc, is therefore 3.64ft -1 co- 
moving Mpc smaller than the difference that would be 
expected under the null hypothesis. A difference as large 
or larger than Aro = 3.64ft -1 Mpc between AGN with 
-25 < M1350 < -19 and -30 < M1350 < -25 occurred 
in 10% of the randomized AGN samples described above. 
We conclude that the null hypothesis of a constant duty 
cycle can be rejected with roughly 90% confidence. 

5. SUMMARY & DISCUSSION 

We measured the galaxy- AGN cross-correlation length 
ro as a function of AGN luminosity. The cross-correlation 
length was similar for bright and faint AGN, ro = 
4.7 ± 2.3 for -30 < Af 135 o < -25 and r = 5.4 ± 1.2 
for —25 < Afi35o < —19, which led us to conclude with 
90% confidence that both are found in halos with simi- 
lar masses and that bright AGN are rarer because their 
duty cycle is shorter. Since halo lifetimes depend only 
weakly on halo mass (e.g., Martini & Weinberg 2001), 
the difference in duty cycle implies that optically faint 
AGN have longer lifetimes. 

Our analysis differs from previous work (e.g., that of 
Croom et al. 2005, who also found no luminosity de- 
pendence in the AGN clustering strength) in two princi- 
pal ways. We estimated the duty cycle from the cross- 
correlation of galaxies and AGN, not from the auto- 
correlation function of AGN, and our sample included 
AGN with a much wider range of luminosities, extend- 
ing ~ 4 magnitudes fainter than the QSO threshold 
Afi35o = —23. These differences allowed us to obtain 
our measurement from a comparatively small survey. 

An appraisal of this result should cover at least the 
following three points. 

The first is obvious: it is only marginally significant. 
Larger samples will be required to prove that the duty 
cycle depends on luminosity. Moreover, other arguments 
suggest that the minimum allowed duty cycle at high lu- 
minosity should be increased and the maximum allowed 
at low luminosity should be decreased. Since the AGN 
lifetime is roughly the age of the universe times the duty 
cycle (e.g., Martini & Weinberg 2001), a duty cycle of 
10 -5 for the brightest AGN is incompatible with the 
observed proximity effect in QSOs' spectra (e.g., Martini 
2004) and with the lack of flickering QSOs in the Sloan 
Digital Sky Survey (Martini & Schneider 2003). A duty 
cycle of roughly unity for the fainter AGN is implausible 
as well, since a black hole radiating continuously would 

6 Note that the lack of a detected AGN in most high-redshift 
galaxies is not by itself an argument against a duty cycle of unity 
for AGN with luminosities -25 < M1350 < -19. These AGN could 



almost certainly be too faint compared to its galaxy 
for us to detect: the difference in energetic efficiency 
for black hole accretion (O.lmc 2 ) and hydrogen burn- 
ing (0.007mc 2 ) implies that a galaxy's steadily radiating 
black hole would be much fainter than its stars if the final 
ratio of black hole to stellar mass is A/bh/A/* ~ 0.001. 
6 Taking these arguments into account would bring the 
high and low luminosity duty cycles closer together in 
Figure El 

Second, the physical interpretation is not straightfor- 
ward. Recall that we have defined the duty cycle for the 
absolute magnitude range M\ Q < M < Mhi as the ratio 
of the number density of AGN with those magnitudes to 
the number density of halos that can host them. In the 
appendix we show that this duty cycle would be indepen- 
dent of magnitude if blackholes accreted only at the Ed- 
dington rate, were not obscured by dust, and had masses 
that followed a tight power-law correlation with the total 
masses Mh of galaxies that contain them. The duty cy- 
cle would decrease at large luminosities if brighter AGN 
were more heavily obscured, if black hole masses fell be- 
low the predictions of the M-QYi-Mh correlation at very 
large Mh, or if anything (e.g., complicated light curves) 
gave a broad range of luminosities L to the AGN that 
lie within halos of a given mass Mh- Each of these is 
expected theoretically (e.g., Hopkins et al. 2005). The 
apparent decrease of duty cycle at large luminosities pre- 
sumably results from a combination of physical effects, 
and our observations do not identify which is dominant 
among them. 

Finally, our result was derived from a small survey de- 
signed for other purposes. Most of the brightest AGN 
lay behind the survey galaxies, not in their midst, reduc- 
ing the number of galaxy- AGN pairs and increasing the 
uncertainty in ro. A large, optimized survey could easily 
shrink the error bars several fold. The only useful con- 
tribution of this paper may be its demonstration that a 
definitive measurement is within easy reach. 
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that they let us proceed with the analysis. This research 
has made use of the NASA/IPAC Extragalactic Database 
(NED) which is operated by the Jet Propulsion Labora- 
tory, California Institute of Technology, under contract 
with the National Aeronautics and Space Administra- 
tion. 

shine exclusively within the most massive galaxies, leaving the less 
massive galaxies with AGN that are undetectably faint. 



APPENDIX 

PHYSICAL INTERPRETATION OF THE DUTY CYCLE 

We discuss three simple models for AGN evolution that may help illustrate the physical meaning of the duty cycle. 
Suppose first that black hole mass Mbh is tightly correlated with total galaxy mass Mh at all times, that the 
correlation has the form Mbh oc Mfi , that AGN are unobscured by dust, and that black holes radiate either at 
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Fig. A7. — Net time spent at a given luminosity for an ensemble of 6 blackholes in the first toy model considered in the appendix. We 
assume that these black holes radiate at the Eddington luminosity and lie inside 6 halos of equal mass M^. If the blackholes have the initial 
luminosities that are marked with vertical lines, and they grow until the luminosities have reached the final values marked with circles, 
then the total amount of time the 6 AGN spent radiating at a given luminosity is shown by the solid histogram. The distribution for all 
AGN in halos of mass M^, not just these 6 AGN, might look more like the dotted curve in the background. If black hole and halo mass 
are tightly correlated and all accretion is at the Eddington rate, this function has to be narrow. The units on the x axis are normalized to 
E(Lp,dd\ Mfr), the mean Eddington luminosity of all AGN in halos of mass units on the y axis are arbitrary. 



the Eddington rate LEdd(MBH) or not at all, gaining their mass in a few short accretion episodes separated by long 
periods of quiescence. The duty cycle would then be independent of AGN magnitude, as can be seen with the following 
argument. 

Begin by considering the evolution of a black hole inside a single dark matter halo of given mass M^. When the halo 
forms in the very early stages of a merger of two smaller halos, its black hole mass 7 may initially be smaller than the 
mean mass implied by the M B n~M h correlation, but by the time the halo is destroyed by mergers, roughly one Hubble 
time later (Martini & Weinberg 2001), the black hole must have grown enough to fall on the correlation. Otherwise 
the correlation could not be satisfied by the ensemble of all halos. Since accretion at the Eddington rate produces 
exponential growth, the black hole will spend equal amounts of time in each octave of luminosity as it grows from 
its initial mass Mi to its final mass M/; if one were to plot the amount of time spent in each logarithmic interval of 
luminosity L, it would be constant for L^ddiMi) < L < LEdd(Mf) and elsewhere. This is equally true if the growth 
occurs in many discrete episodes of accretion or in a single burst. Now consider a plot of total elapsed time versus 
luminosity for the black holes within N randomly chosen halos of the same mass M^. It would be the superposition of 
N boxcars with random left and right edges, producing an overall shape that is peaked near the Eddington luminosity 
of the typical black hole associated with halos of mass Mh- Figure lATI shows an example for N = 6. The same plot 
for the ensemble of all halos of mass Mh would be a smoother realization of a similar function. Call this plot the 
kernel. Since the number of AGN we observe with a given luminosity is proportional to the net time AGN spend at 
that luminosity, the kernel is the AGN luminosity distribution we would observe if the universe consisted solely of 
halos with mass Mh- The width of the kernel depends on how far the initial and final black hole masses stray from 
the expectation value E(MBn\Mh), but it must be very narrow compared to the multi-decade width of the halo mass 
distribution. Otherwise our assumption of a tight Mbh - -Mfc correlation would be violated. The AGN within a narrow 
range of luminosity therefore must lie inside halos with a narrow range of masses. Our definition of duty cycle for 
L m in < L < i max is the number density of AGN within that range of luminosity divided by the number density of 
halos that can host them. In this scenario, it is equal to the time required for the AGN's luminosity to grow from I/ m i n 
to £ ma x if it is accreting at the Eddington rate divided by the halos' mean lifetime. The numerator is independent of 
halo mass for logarithmic luminosity intervals, and the denominator depends extremely weakly on halo mass (Martini 
& Weinberg 2001). Therefore the duty cycle in logarithmic luminosity bins should be nearly independent of halo mass 
or AGN luminosity. 

To check this claim, we generated an ensemble of simulated AGN by starting with an ensemble of halos following 
a Press-Schechter mass function (CIm = 0.3, = 0.7, T = 0.2, erg = 0.9, z = 2.5), assigning each halo an expected 

7 which may initially be divided among two black holes; since the Eddington luminosity of black hole of mass 2M is equal to the sum of 
the Eddington luminosities two black holes each of mass M, this does not affect our argument 
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central blackhole mass with the relationship M BH = 10 7 (M /l /10 12 M Q ) 1 - 65 (Ferrarese 2002), and giving each black 
hole a luminosity equal t o th e Eddington luminosity of the expected mass times a random number drawn at random 
from t he k ernel (Figure IA7|I . This resulted in the AGN luminosity distribution shown in the upper left panel of 
Figure IA8I The distribution of halo masses for AGN with luminosities L E dd(10 6 Af o ) < L < L E dd(lO 8 M ) and 
iEdd(lO 8 M0) < L < LEdd(10 10 MQ) is shown in the middle left panel. The bottom panel shows the inferred duty 
cycle in these luminosity ranges, i.e., the ratio of AGN number density in each luminosity range to the number density 
of halos more massive than the mean associated halo mass shown in the middle left panel. This is roughly the duty 
cycle that would be estimated with the approach we adopted above. It is the same for the two logarithmic luminosity 
ranges, as expected. 

The scenario can be altered in two ways to make the duty cycle decrease at larger luminosities. 

The first is to increase the number density of halos that can host the brightest AGN. For a fixed halo mass distribution, 
this can be accomplished by relaxing our assumption that accretion occurs only at the Eddington rate or by increasing 
the scatter in the Msa-Mh relationship. Either increases the scatter in the relationship between Mh and L, r aising 
the probability that a high luminosity AGN resides within a low mass halo. The middle column of Figure lA"8l shows 
one example of how a broad distribution of L at fixed Mh makes the duty cycle depend on luminosity. 

The second way is to reduce the lifetimes of the brightest AGN. If the M^K-Mh correlation is a tight power-law and all 
accretion is at the Eddington rate, then we can adjust neither the mean accreted mass for blackholes in the most massive 
halos nor the rate at which accretion occurs. In this case lifetimes of the brightest AGN can be reduced only by making 
them heavily obscured while they accrete most of their mass. The lifetimes can also be reduced, even for unobscured 
Eddington-rate accretion, if we change the form of the M^n-Mh relationship. One change seems well motivated: letting 
it break down for halos with super-galactic masses. Ferrarese's (2002) relationship M BH ~ 10 7 (M /l /10 12 M Q ) 1 - 65 M Q 
predicts that local clusters of mass 10 15 M Q should contain 10 12 M Q central blackholes, for example, but there is no 
evidence that these ultra-massive blackhole exist. It seems more likely that blackhole formation becomes as suppressed 
as star- formation in halos with mass Mh 10 13 M Q . Suppressing or obscuring the brightest AGN can make the duty 
cycle depend strongly on luminosity, as Figure lA~8l shows . 

If additional observations confirm the decrease in duty cycle at high luminosities, some combination of these effects 
would presumably be responsible. 
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Fig. A8. — Dependence of duty cycle on luminosity for three toy models. We generated an ensemble of simulated central black holes from 
a Prcss-Schcchtcr mass function by associating each halo of mass with a blackhole of expected mass Mbh/Mq = 10 7 (Mj 1 /10 12 Mq) 1 ' 65 
(Ferrarese 2002), then associated each blackhole with a luminosity under different assumptions for the three models. Each column shows 
results for one model. Left panels assume Eddington accretion and a tight correlation of and Mbh- Middle panels assume large scatter 
in the relationship between M h and AGN luminosity L. Right panels assume Eddington accretion but stunt the growth of blackholes in 
the most massive halos. Units on the y axis are arbitrary in all panels. Top panels: Dashed lines show the luminosity distribution that 
would have resulted if each black hole had exactly the mean mass predicted by the Mrh-M^ relationship and radiated at its Eddington 
luminosity. Instead we assumed that each halo's AGN luminosity had some scatter around its expectation value. Dotted lines indicate 
the assumed luminosity distribution for halos of mass 10 12 , 10 13 , 10 14 , 10 15 Mq. They show the assumed scatter in the M^-L relationship, 
which is large for the middle model and equally small for the left and right models. The solid lines show the implied AGN luminosity 
function. (More-realistic models would keep the AGN luminosity function fixed to the observations by adjusting other parameters, but 
these simple models are sufficient to illustrate our point.) Middle panels: Distribution of mass for halos that host AGN with luminosities 
10 6 L|, dd < L < 10 8 L| dd and 10 8 L| dd < L < 10 lo L| dd . Here L| dd is the Eddington luminosity of a solar-mass blackhole. Bottom 
panels: Derived duty cycle for AGN in the same luminosity ranges. The duty cycles were estimated by dividing the number of AGN in 
the luminosity range by the number density of halos more massive than the mean shown in the middle panels. These panels show that the 
duty cycle will decrease at high luminosities if there is significant scatter in the M^-L relationship or if black hole accretion is suppressed 
in the most massive halos. Increases in dust obscuration with luminosity can also reduce the duty cycle at high L, but wc judged this effect 
too obvious to illustrate here. 



